Variant mutation in SARS-CoV-2 nucleocapsid enhances viral infection via altered genomic encapsidation (preprint)

The evolution of SARS-CoV-2 variants and their respective phenotypes represents an important set of tools to understand basic coronavirus biology as well as the public health implications of individual mutations in variants of concern. While mutations outside of Spike are not well studied, the entire viral genome is undergoing evolutionary selection, particularly the central disordered linker region of the nucleocapsid (N) protein. Here, we identify a mutation (G215C), characteristic of the Delta variant, that introduces a novel cysteine into this linker domain, which results in the formation of a disulfide bond and a stable N-N dimer. Using reverse genetics, we determined that this cysteine residue is necessary and sufficient for stable dimer formation in a WA1 SARS-CoV-2 background, where it results in significantly increased viral growth both in vitro and in vivo. Finally, we demonstrate that the N:G215C virus packages more nucleocapsid per virion and that individual virions are larger, with elongated morphologies.

SARS-CoV-2 hijacks fragile X mental retardation proteins for efficient infection (preprint)

Viruses interact with numerous host factors to facilitate viral replication and to dampen antiviral defense mechanisms. We currently have a limited mechanistic understanding of how SARS-CoV-2 binds host factors and the functional role of these interactions. Here, we uncover a novel interaction between the viral NSP3 protein and the fragile X mental retardation proteins (FMRPs: FMR1 and FXR1-2). SARS-CoV-2 NSP3 mutant viruses preventing FMRP binding have attenuated replication in vitro and have delayed disease onset in vivo. We show that a unique peptide motif in NSP3 binds directly to the two central KH domains of FMRPs and that this interaction is disrupted by the I304N mutation found in a patient with fragile X syndrome. NSP3 binding to FMRPs disrupts their interaction with the stress granule component UBAP2L through direct competition with a peptide motif in UBAP2L to prevent FMRP incorporation into stress granules. Collectively, our results provide novel insight into how SARS-CoV-2 hijacks host cell proteins for efficient infection and provides molecular insight to the possible underlying molecular defects in fragile X syndrome.

Interaction between host G3BP and viral nucleocapsid protein regulates SARS-CoV-2 replication (preprint)

G3BP1/2 are paralogous proteins that promote stress granule formation in response to cellular stresses, including viral infection. G3BP1/2 are prominent interactors of the nucleocapsid (N) protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). However, the functional consequences of the G3BP1-N interaction in the context of viral infection remain unclear. Here we used structural and biochemical analyses to define the residues required for G3BP1-N interaction, followed by structure-guided mutagenesis of G3BP1 and N to selectively and reciprocally disrupt their interaction. We found that mutation of F17 within the N protein led to selective loss of interaction with G3BP1 and consequent failure of the N protein to disrupt stress granule assembly. Introduction of SARS-CoV-2 bearing an F17A mutation resulted in a significant decrease in viral replication and pathogenesis in vivo, indicating that the G3BP1-N interaction promotes infection by suppressing the ability of G3BP1 to form stress granules.

The impact of RSV/SARS-CoV-2 co-infection on clinical disease and viral replication: insights from a BALB/c mouse model (preprint)

RSV and SARS-CoV-2 are prone to co-infection with other respiratory viruses. In this study, we use RSV/SARS-CoV-2 co-infection to evaluate changes to clinical disease and viral replication in vivo. To consider the severity of RSV infection, effect of sequential infection, and the impact of infection timing, mice were co-infected with varying doses and timing. Compared with a single infection of RSV or SARS-CoV-2, the co-infection of RSV/SARS-CoV-2 and the primary infection of RSV followed by SARS-CoV-2 results in protection from SARS-CoV-2-induced clinical disease and reduces SARS-CoV-2 replication. Co-infection also augmented RSV replication at early timepoints with only the low dose. Additionally, the sequential infection of RSV followed by SARS-CoV-2 led to improved RSV clearance regardless of viral load. However, SARS-CoV-2 infection followed by RSV results in enhanced SARS-CoV-2-induced disease while protecting from RSV-induced disease. SARS-CoV-2/RSV sequential infection also reduced RSV replication in the lung tissue, regardless of viral load. Collectively, these data suggest that RSV and SARS-CoV-2 co-infection may afford protection from or enhancement of disease based on variation in infection timing, viral infection order, and/or viral dose. In the pediatric population, understanding these infection dynamics will be critical to treat patients and mitigate disease outcomes.

Loss-of-function mutation in Omicron variants reduces spike protein expression and attenuates SARS-CoV-2 infection (preprint)

SARS-CoV-2 Omicron variants emerged in 2022 with >30 novel amino acid mutations in the spike protein alone. While most studies focus on the impact of receptor binding domain changes, mutations in the C-terminal of S1 (CTS1), adjacent to the furin cleavage site, have largely been ignored. In this study, we examined three Omicron mutations in CTS1: H655Y, N679K, and P681H. Generating a SARS-CoV-2 triple mutant (YKH), we found that the mutant increased spike processing, consistent with prior reports for H655Y and P681H individually. Next, we generated a single N679K mutant, finding reduced viral replication in vitro and less disease in vivo. Mechanistically, the N679K mutant had reduced spike protein in purified virions compared to wild-type; spike protein decreases were further exacerbated in infected cell lysates. Importantly, exogenous spike expression also revealed that N679K reduced overall spike protein yield independent of infection. Together, the data show that N679K is a loss-of-function mutation reducing overall spike levels during omicron infection, which may have important implications for disease severity, immunity, and vaccine efficacy.

SARS-CoV-2 Uses Nonstructural Protein 16 to Evade Restriction by IFIT1 and IFIT3 (preprint)

Understanding the molecular basis of innate immune evasion by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an important consideration for designing the next wave of therapeutics. Here, we investigate the role of the nonstructural protein 16 (NSP16) of SARS-CoV-2 in infection and pathogenesis. NSP16, a ribonucleoside 2'-O methyltransferase (MTase), catalyzes the transfer of a methyl group to mRNA as part of the capping process. Based on observations with other CoVs, we hypothesized that NSP16 2'-O MTase function protects SARS-CoV-2 from cap-sensing host restriction. Therefore, we engineered SARS-CoV-2 with a mutation that disrupts a conserved residue in the active site of NSP16. We subsequently show that this mutant is attenuated both in vitro and in vivo, using a hamster model of SARS-CoV-2 infection. Mechanistically, we confirm that the NSP16 mutant is more sensitive to type I interferon (IFN-I) in vitro. Furthermore, silencing IFIT1 or IFIT3, IFN-stimulated genes that sense a lack of 2'-O methylation, partially restores fitness to the NSP16 mutant. Finally, we demonstrate that sinefungin, a methyltransferase inhibitor that binds the catalytic site of NSP16, sensitizes wild-type SARS-CoV-2 to IFN-I treatment. Overall, our findings highlight the importance of SARS-CoV-2 NSP16 in evading host innate immunity and suggest a possible target for future antiviral therapies.

The Arabidopsis Framework Model version 2 predicts the organism-level effects of circadian clock gene mis-regulation (preprint)

Predicting a multicellular organism's phenotype quantitatively from its genotype is challenging, as genetic effects must propagate across scales. Circadian clocks are intracellular regulators that control temporal gene expression patterns and hence metabolism, physiology and behaviour. Here we explain and predict canonical phenotypes of circadian timing in a multicellular, model organism. We used diverse metabolic and physiological data to combine and extend mathematical models of rhythmic gene expression, photoperiod-dependent flowering, elongation growth and starch metabolism within a Framework Model for the vegetative growth of Arabidopsis thaliana, sharing the model and data files in a structured, public resource. The calibrated model predicted the effect of altered circadian timing upon each particular phenotype in clock-mutant plants under standard laboratory conditions. Altered night-time metabolism of stored starch accounted for most of the decrease in whole-plant biomass, as previously proposed. Mobilisation of a secondary store of malate and fumarate was also mis-regulated, accounting for any remaining biomass defect. We test three candidate mechanisms for the accumulation of these organic acids. Our results link genotype through specific processes to higher-level phenotypes, formalising our understanding of a subtle, pleiotropic syndrome at the whole-organism level, and validating the systems approach to understand complex traits starting from intracellular circuits. This work updates the first biorXiv version, February 2017, https://doi.org/10.1101/105437, with an expanded description and additional analysis of the same core data sets and the same FMv2 model, summary tables and supporting, follow-on data from three further studies with further collaborators. This biorXiv revision constitutes the second version of this report.

A live-attenuated SARS-CoV-2 vaccine candidate with accessory protein deletions (preprint)

We report a live-attenuated SARS-CoV-2 vaccine candidate with (i) re-engineered viral transcriptional regulator sequences and (ii) deleted open-reading-frames (ORF) 3, 6, 7, and 8 ({triangleup}3678). The {triangleup}3678 virus replicates about 7,500-fold lower than wild-type SARS-CoV-2 on primary human airway cultures, but restores its replication on interferon-deficient Vero-E6 cells that are approved for vaccine production. The {triangleup}3678 virus is highly attenuated in both hamster and K18-hACE2 mouse models. A single-dose immunization of the {triangleup}3678 virus protects hamsters from wild-type virus challenge and transmission. Among the deleted ORFs in the {triangleup}3678 virus, ORF3a accounts for the most attenuation through antagonizing STAT1 phosphorylation during type-I interferon signaling. We also developed an mNeonGreen reporter {triangleup}3678 virus for high-throughput neutralization and antiviral testing. Altogether, the results suggest that {triangleup}3678 SARS-CoV-2 may serve as a live-attenuated vaccine candidate and a research tool for potential biosafety level-2 use.

Engineering SARS-CoV-2 cocktail antibodies into a bispecific format improves neutralizing potency and breadth (preprint)

One major limitation of neutralizing antibody-based COVID-19 therapy is the requirement of costly cocktails to reduce antibody resistance. We engineered two bispecific antibodies (bsAbs) using distinct designs and compared them with parental antibodies and their cocktail. Single molecules of both bsAbs block the two epitopes targeted by parental antibodies on the receptor-binding domain (RBD). However, bsAb with the IgG-(scFv)2 design (14-H-06) but not the CrossMAb design (14-crs-06) increases antigen-binding and virus-neutralizing activities and spectrum against multiple SARS-CoV-2 variants including the Omicron, than the cocktail. X-ray crystallography and computational simulations reveal distinct neutralizing mechanisms for individual cocktail antibodies and suggest higher inter-spike crosslinking potentials by 14-H-06 than 14-crs-06. In mouse models of infections by SARS-CoV-2 and the Beta, Gamma, and Delta variants, 14-H-06 exhibits higher or equivalent therapeutic efficacy than the cocktail. Rationally engineered bsAbs represent a cost-effective alternative to antibody cocktails and a promising strategy to improve potency and breadth.

QTQTN motif upstream of the furin-cleavage site plays key role in SARS-CoV-2 infection and pathogenesis. (preprint)

The furin cleavage site (FCS), an unusual feature in the SARS-CoV-2 spike protein, has been spotlighted as a factor key to facilitating infection and pathogenesis by increasing spike processing 1,2. Similarly, the QTQTN motif directly upstream of the FCS is also an unusual feature for group 2B coronaviruses (CoVs). The QTQTN deletion has consistently been observed in in vitro cultured virus stocks and some clinical isolates 3. To determine whether the QTQTN motif is critical to SARS-CoV-2 replication and pathogenesis, we generated a mutant deleting the QTQTN motif ({Delta}QTQTN). Here we report that the QTQTN deletion attenuates viral replication in respiratory cells in vitro and attenuates disease in vivo. The deletion results in a shortened, more rigid peptide loop that contains the FCS, and is less accessible to host proteases, such as TMPRSS2. Thus, the deletion reduced the efficiency of spike processing and attenuates SARS-CoV-2 infection. Importantly, the QTQTN motif also contains residues that are glycosylated4, and disruption its glycosylation also attenuates virus replication in a TMPRSS2-dependent manner. Together, our results reveal that three aspects of the S1/S2 cleavage site (the FCS, loop length, and glycosylation) are required for efficient SARS-CoV-2 replication and pathogenesis.

Nucleocapsid mutations in SARS-CoV-2 augment replication and pathogenesis. (preprint)

While SARS-CoV-2 continues to adapt for human infection and transmission, genetic variation outside of the spike gene remains largely unexplored. This study investigates a highly variable region at residues 203-205 in SARS-CoV-2 nucleocapsid protein. Recreating the alpha variant mutation in an early pandemic (WA-1) background, we found that the R203K-G204R mutation is sufficient to enhance replication, fitness, and pathogenesis of SARS-CoV-2. Importantly, the R203K-G204R mutation increases nucleocapsid phosphorylation, providing a molecular basis for these phenotypes. Notably, an analogous alanine substitution mutant also increases SARS-CoV-2 fitness and phosphorylation, suggesting that infection is enhanced through ablation of the ancestral RG motif. Overall, these results demonstrate that variant mutations outside spike are also key components in SARS-CoV-2 continued adaptation to human infection.

Delta spike P681R mutation enhances SARS-CoV-2 fitness over Alpha variant (preprint)

SARS-CoV-2 Delta variant has rapidly replaced the Alpha variant around the world. The mechanism that drives this global replacement has not been defined. Here we report that Delta spike mutation P681R plays a key role in the Alpha-to-Delta variant replacement. In a replication competition assay, Delta SARS-CoV-2 efficiently outcompeted the Alpha variant in human lung epithelial cells and primary human airway tissues. Delta SARS-CoV-2 bearing the Alpha-spike glycoprotein replicated less efficiently than the wild-type Delta variant, suggesting the importance of Delta spike in enhancing viral replication. The Delta spike has accumulated mutation P681R located at a furin cleavage site that separates the spike 1 (S1) and S2 subunits. Reverting the P681R mutation to wild-type P681 significantly reduced the replication of Delta variant, to a level lower than the Alpha variant. Mechanistically, the Delta P681R mutation enhanced the cleavage of the full-length spike to S1 and S2, leading to increased infection via cell surface entry. In contrast, the Alpha spike also has a mutation at the same amino acid (P681H), but the spike cleavage from purified Alpha virions was reduced compared to the Delta spike. Collectively, our results indicate P681R as a key mutation in enhancing Delta variant replication via increased S1/S2 cleavage. Spike mutations that potentially affect furin cleavage efficiency must be closely monitored for future variant surveillance.

Common Mechanism of SARS-CoV and SARS-CoV-2 Pathogenesis across Species (preprint)

Sarbecovirus (CoV) infections, including Severe Acute Respiratory CoV (SARS-CoV) and SARS-CoV-2, are considerable human threats. Human GWAS studies have recently identified loci associated with variation in SARS-CoV-2 susceptibility. However, genetically tractable models that reproduce human CoV disease outcomes are needed to mechanistically evaluate genetic determinants of CoV susceptibility. We used the Collaborative Cross (CC) and human GWAS datasets to elucidate host susceptibility loci that regulate CoV infections and to identify host quantitative trait loci that modulate severe CoV and pan-CoV disease outcomes including a major disease regulating loci including CCR9. CCR9 ablation resulted in enhanced titer, weight loss, respiratory dysfunction, mortality, and inflammation, providing mechanistic support in mitigating protection from severe SARS-CoV-2 pathogenesis across species. This study represents a comprehensive analysis of susceptibility loci for an entire genus of human pathogens conducted, identifies a large collection of susceptibility loci and candidate genes that regulate multiple aspects type-specific and cross-CoV pathogenesis, and also validates the paradigm of using the CC platform to identify common cross-species susceptibility loci and genes for newly emerging and pre-epidemic viruses.

Mouse Adapted SARS-CoV-2 protects animals from lethal SARS-CoV challenge (preprint)

The emergence of SARS-CoV-2 has resulted in a worldwide pandemic causing significant damage to public health and the economy. Efforts to understand the mechanisms of COVID-19 disease have been hampered by the lack of robust mouse models. To overcome this barrier, we utilized a reverse genetic system to generate a mouse-adapted strain of SARS-CoV-2. Incorporating key mutations found in SARSCoV-2 variants, this model recapitulates critical elements of human infection including viral replication in the lung, immune cell infiltration, and significant in vivo disease. Importantly, mouse-adaptation of SARS-CoV-2 does not impair replication in human airway cells and maintains antigenicity similar to human SARS-CoV-2 strains. Utilizing this model, we demonstrate that SARS-CoV-2 infected mice are protected from lethal challenge with the original SARS-CoV, suggesting immunity from heterologous CoV strains. Together, the results highlight the utility of this mouse model for further study of SARS CoV-2 infection and disease.

CCR2-dependent monocyte-derived cells restrict SARS-CoV-2 infection (preprint)

SARS-CoV-2 has caused a historic pandemic of respiratory disease (COVID-19) and current evidence suggests severe disease is associated with dysregulated immunity within the respiratory tract. However, the innate immune mechanisms that mediate protection during COVID-19 are not well defined. Here we characterize a mouse model of SARS-CoV-2 infection and find that early CCR2-dependent infiltration of monocytes restricts viral burden in the lung. We find that a recently developed mouse-adapted MA-SARS-CoV-2 strain, as well as the emerging B.1.351 variant, trigger an inflammatory response in the lung characterized by expression of pro-inflammatory cytokines and interferon-stimulated genes. scRNA-seq analysis of lung homogenates identified a hyper-inflammatory monocyte profile. Using intravital antibody labeling, we demonstrate that MA-SARS-CoV-2 infection leads to increases in circulating monocytes and an influx of CD45+ cells into the lung parenchyma that is dominated by monocyte-derived cells. We utilize this model to demonstrate that mechanistically, CCR2 signaling promotes infiltration of classical monocytes into the lung and expansion of monocyte-derived cells. Parenchymal monocyte-derived cells appear to play a protective role against MA-SARS-CoV-2, as mice lacking CCR2 showed higher viral loads in the lungs, increased lung viral dissemination, and elevated inflammatory cytokine responses. These studies have identified a CCR2-monocyte axis that is critical for promoting viral control and restricting inflammation within the respiratory tract during SARS-CoV-2 infection.

Nucleocapsid vaccine elicits spike-independent SARS-CoV-2 protective immunity (preprint)

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the COVID-19 pandemic. Neutralizing antibodies target the receptor binding domain of the spike (S) protein, a focus of successful vaccine efforts. Concerns have arisen that S-specific vaccine immunity may fail to neutralize emerging variants. We show that vaccination with HAd5 expressing the nucleocapsid (N) protein can establish protective immunity, defined by reduced weight loss and viral load, in both Syrian hamsters and k18-hACE2 mice. Challenge of vaccinated mice was associated with rapid N-specific T cell recall responses in the respiratory mucosa. This study supports the rationale for including additional viral antigens, even if they are not a target of neutralizing antibodies, to broaden epitope coverage and immune effector mechanisms.

Endomembrane targeting of human OAS1 p46 augments antiviral activity (preprint)

Many host RNA sensors are positioned in the cytosol to detect viral RNA during infection. However, most positive-strand RNA viruses replicate within a modified organelle co-opted from intracellular membranes of the endomembrane system, which shields viral products from host cell innate immune sensors. Targeting innate RNA sensors to the endomembrane system may enhance their ability to sense viral RNA generated by viruses that use these compartments for replication. Here, we reveal that an isoform of oligoadenylate synthetase 1, OAS1 p46, is prenylated and targeted to the endomembrane system. Membrane localization of OAS1 p46 confers enhanced access to viral replication sites and results in increased antiviral activity against a subset of RNA viruses including flavivirus, picornavirus, and SARS-CoV-2. Finally, our human genetic analysis shows that the OAS1 splice-site SNP responsible for production of the OAS1 p46 isoform strongly associates with COVID-19 severity. This study highlights the importance of endomembrane targeting for the antiviral specificity of OAS1 and suggests early control of SARS-CoV-2 replication through OAS1-p46 is an important determinant of COVID-19 severity.

Tiled-ClickSeq for targeted sequencing of complete coronavirus genomes with simultaneous capture of RNA recombination and minority variants (preprint)

High-throughput genomics of SARS-CoV-2 is essential to characterize virus evolution and to identify adaptations that affect pathogenicity or transmission. While single-nucleotide variations (SNVs) are commonly considered as driving virus adaption, RNA recombination events that delete or insert nucleic acid sequences are also critical. Whole genome targeting sequencing of SARS-CoV-2 is typically achieved using pairs of primers to generate cDNA amplicons suitable for Next-Generation Sequencing (NGS). However, paired-primer approaches impose constraints on where primers can be designed, how many amplicons are synthesized and requires multiple PCR reactions with non-overlapping primer pools. This imparts sensitivity to underlying SNVs and fails to resolve RNA recombination junctions that are not flanked by primer pairs. To address these limitations, we have designed an approach called 'Tiled-ClickSeq'. Tiled-ClickSeq uses hundreds of tiled-primers spaced evenly along the virus genome in a single reverse-transcription reaction. The other end of the cDNA amplicon is generated by azido-nucleotides that stochastically terminate cDNA synthesis, obviating the need for a paired-primer. A sequencing adaptor containing a Unique Molecular Identifier (UMI) is appended using click-chemistry and a PCR reaction using Illumina adaptors generates a final NGS library. Tiled-ClickSeq provides complete genome coverage, including the 5'UTR, at high depth and specificity to virus on both Illumina and Nanopore NGS platforms. Here, we analyze multiple SARS-CoV-2 isolates and simultaneously characterize minority variants, sub-genomic mRNAs (sgmRNAs), structural variants (SVs) and D-RNAs. Tiled-ClickSeq therefore provides a convenient and robust platform for SARS-CoV-2 genomics that captures the full range of RNA species in a single, simple assay.

The N501Y spike substitution enhances SARS-CoV-2 transmission (preprint)

Beginning in the summer of 2020, a variant of SARS-CoV-2, the cause of the COVID-19 pandemic, emerged in the United Kingdom (UK). This B.1.1.7 variant increased rapidly in prevalence among sequenced strains, attributed to an increase in infection and/or transmission efficiency. The UK variant has 19 nonsynonymous mutations across its viral genome including 8 substitutions or deletions in the spike protein, which interacts with cellular receptors to mediate infection and tropism. Here, using a reverse genetics approach, we show that, of the 8 individual spike protein substitutions, only N501Y exhibited consistent fitness gains for replication in the upper airway in the hamster model as well as primary human airway epithelial cells. The N501Y substitution recapitulated the phenotype of enhanced viral transmission seen with the combined 8 UK spike mutations, suggesting it is a major determinant responsible for increased transmission of this variant. Mechanistically, the N501Y substitution improved the affinity of the viral spike protein for cellular receptors. As suggested by its convergent evolution in Brazil and South Africa, our results indicate that N501Y substitution is a major adaptive spike mutation of major concern.

Neutralization of spike 69/70 deletion, E484K, and N501Y SARS-CoV-2 by BNT162b2 vaccine-elicited sera (preprint)

We engineered three SARS-CoV-2 viruses containing key spike mutations from the newly emerged United Kingdom (UK) and South African (SA) variants: N501Y from UK and SA; 69/70-deletion+N501Y+D614G from UK; and E484K+N501Y+D614G from SA. Neutralization geometric mean titers (GMTs) of twenty BTN162b2 vaccine-elicited human sera against the three mutant viruses were 0.81- to 1.46-fold of the GMTs against parental virus, indicating small effects of these mutations on neutralization by sera elicited by two BNT162b2 doses.